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In a Bose superfluid, the coupling between transverse (phase) and longitudinal fluctuations leads 
to a divergence of the longitudinal correlation function, which is responsible for the occurrence of 
infrared divergences in the perturbation theory and the breakdown of the Bogoliubov approximation. 
We report a non-perturbative renormalization-group (NPRG) calculation of the one-particle Green 
function of an interacting boson system at zero temperature. We find two regimes separated by a 
characteristic momentum scale ka ("Ginzburg" scale). While the Bogoliubov approximation is valid 
at large momenta and energies, |p|, \u)\/c 3> ka (with c the velocity of the Bogoliubov sound mode), 
in the infrared (hydrodynamic) regime |p|,|cj|/c <C ka the normal and anomalous self-energies 
exhibit singularities reflecting the divergence of the longitudinal correlation function. In particular, 
we find that the anomalous self-energy agrees with the Bogoliubov result E an (p,cj) ~ const at high- 
energies and behaves as E an (p,o;) ~ (c 2 p 2 — uj 2 )^ d ~ 3 ^ 2 in the infrared regime (with d the space 
dimension), in agreement with the Nepomnyashchii identity E an (0, 0) = and the predictions of 
Popov's hydrodynamic theory. We argue that the hydrodynamic limit of the one-particle Green 
function is fully determined by the knowledge of the exponent 3 — d characterizing the divergence of 
the longitudinal susceptibility and the Ward identities associated to gauge and Galilean invariances. 
The infrared singularity of E an (p,u) leads to a continuum of excitations (coexisting with the sound 
mode) which shows up in the one-particle spectral function. 

PACS numbers: 05.30.Jp,03.75.Kk,05.10.Cc 



I. INTRODUCTION 

Following the success of the Bogoliubov theory [l[ 
in providing a microscopic explanation of superfluidity, 
much theoretical work has been devoted to the under- 
standing of the infrared behavior and the calculation of 
the one-particle Green function of a Bose superfluid Q. 
Early attempts to improve the Bogoliubov approxima- 
tion encountered difficulties due to a singularperturba- 
tion theory plagued by infrared divergences [|[ 0, H, Hi • 
Although these divergences cancel out in local gauge 
invariant physical quantities (condensate density, Gold- 
stone mode velocity, etc.), they do have a definite physi- 
cal origin: they reflect the divergence of the longitudinal 
susceptibility which is induced by the coupling between 
transverse (phase) and longitudinal fluctuations. This is 
a general phenomenon 0] in systems with a continuous 
broken symmetry as discussed at the end of this section. 

Using field-theoretical diagrammatic methods to han- 
dle the infrared divergences of the perturbation theory, 
Nepomnyashchii and Nepomnyashchii (NN) showed that 
one of the fundamental quantities of a Bose superfluid, 
the anomalous self-energy £ an (p), vanishes in the limit 
p = (p, lu) — ► in dimension d < 3, even though the 
low-energy mode remains phonon-like with a linear spec- 
trum d,H, 03 ■ This exact result shows that the Bo- 
goliubov approximation, where the linear spectrum and 
the superfluidity rely on a finite value of the anomalous 
self-energy, breaks down at low energy. 

An alternative approach to superfluidity, based on 
a phase-amplitude representation of the boson field, 
has been proposed by Popov fill ] . This approach is 



free of infrared singularity, but restricted to the (low- 
momentum) hydrodynamic regime and therefore does 
not allow to study the high- momentum or high-frequency 
regime where the Bogoliubov approximation is expected 
to be valid. Nevertheless, Popov's theory [l2j agrees with 
the asymptotic low-energy behavior of £ an (p) obtained 
by NN 0, 0, Furthermore, the expression of the 

anomalous self-energy obtained by NN and Popov in the 
low-energy limit yields a continuum of (one-particle) exci- 
tations coexisting with the Bogoliubov sound mode [l3T ] , 
in marked contrast with the Bogoliubov theory where the 
sound mode is the sole excitation at low energy. 

The instability of the Bogoliubov fixed point in dimen- 
sion d < 3 towards a different fixed point characterized by 
the divergence of the longitudinal susceptibility has been 
confirmed by Castellani et al. [14|, [l5| . Using the Ward 
identities associated to the local gauge symmetry and a 
renormalization group approach, these authors obtained 
the exact infrared behavior of a Bose superfluid at zero 
temperature. Related results, both at zero flU [171. [l8l. [l9j 
and finite [2(| Hi], HH, [H, [24[ temperature have been ob- 
tained by several authors within the framework of the 
non-perturbative renormalization group. 

In this paper, we study a weakly interacting Bose su- 
perfluid at zero temperature using the so-called BMW 
NPRG scheme introduced by Blaizot, Mcndez-Galain 
and Wschebor [25|, [26|. Compared to more traditional 
RG approaches, the NPRG approach presents a number 
of advantages: i) symmetries are naturally implemented 
(by a proper Ansatz for the effective action or the two- 
point vertex (Sec. |nj)) and Ward identities are naturally 
satisfied; ii) the NPRG approach is not restricted to the 



2 



low-energy asymptotic behavior but can deal with all en- 
ergy scales. In particular, it relates physical quantities at 
a macroscopic scale to the parameters of the microscopic 
model; iii) the BMW scheme enables to obtain the full 
momentum and frequency dependence of the correlation 
functions [l9j |. 

In the NPRG approach, the main quantities of interest 
are the average effective action Tk (the generating func- 
tional of one-particle irreducible vertices) and its SeCOnd- 
^n^ 

order functional derivative, the two-point vertex T k (p), 
whose inverse gives the one-particle propagator (Sec. HI)) . 
Fluctuations beyond the Bogoliubov approximation are 
gradually taken into account as the (running) momentum 
scale k is lowered from the microscopic scale A down to 
zero. In Sec. IIII1 we show that the infrared behavior 
of the one-particle propagator is entirely determined by 
the Ward identities associated to the Galilean and lo- 
cal gauge invariances of the microscopic action, and the 
exponent 3 — d characterizing the divergence of the longi- 
tudinal susceptibility (see the discussion at the end of the 
Introduction). In Sec. HVI wc derive the BMW flow equa- 
tions satisfied by the two-point vertex (p) and obtain 
the analytical solution in the infrared regime. Numer- 
ical results are discussed in Sec. [V] We find that the 
Bogoliubov approximation breaks down at a characteris- 
tic momentum length k(j ("Ginzburg" scale) which, for 
weak boson-boson interactions, is much smaller than the 
inverse healing length kh (kh is defined in Appendix fA"| . 
Although local gauge invariant quantities are not sensi- 
tive to fcc, the effective action Tk is attracted to a fixed 
point characterized by the divergence of the longitudinal 
susceptibility when k kc- We discuss in detail the 
frequency and momentum dependence of the two-point 
vertex T { k 2) (p). While for |p| > k G or \uj\/c > k G (with 

c the velocity of the Bogoliubov sound mode), T k 2 ^(p) is 
well described by the Bogoliubov approximation, we re- 
produce the low-energy asymptotic behavior obtained by 
NN when |p|, \cj\/c -C &g- In this regime, the longitudi- 
nal correlation function becomes singular and its spectral 
function exhibits a continuum of one-particle excitations 
in agreement with the predictions of Popov's hydrody- 
namic theory. Thus our approach provides a unified pic- 
ture of superfluidity in interacting boson systems and 
connects Bogoliubov's theory to Popov's hydrodynamic 
theory. In the conclusion, we comment about a possible 
extension of our results to strongly interacting or one- 
dimensional supcrfluids. 



Infrared behavior of interacting bosons 

Since the divergence of the longitudinal susceptibility 
plays a key role in the infrared behavior of interacting 
boson systems, we first discuss this phenomenon both in 
classical and quantum systems (for a pedagogical discus- 
sion, see also Ref. [27j ). Let us first consider a tp 4 theory 



defined by the action 

s = \l ^{(v^) 2 + V + ^ 2 ) 2 }, (i) 

where cp(r) is a real A^-component field and d the space 
dimension. When v < 0, the mean-field (saddle-point) 
analysis predicts a non-zero order parameter ip = (ip). 
Including Gaussian fluctuations about the saddle-point 
<p Q , we find a gapped mode and N — 1 Goldstonc modes 
corresponding to longitudinal (Sip \\ tp ) and transverse 
(Sip _L ip Q ) fluctuations, respectively The correlation 
functions read 

p +M, (2) 
G±(p) = pV 

This result, which neglects interactions between longitu- 
dinal and transverse fluctuations, is incorrect. In the or- 
dered phase, the amplitude fluctuations of <p are gapped 
and the low-energy effective description is a non-linear a 
model H| 

S[n] = IJ d d r(Vn)\ (3) 

where n is a unit vector (n 2 = 1). To a first approxima- 
tion, equation ([3]) can be obtained by setting ip = |¥>ol n 
in |T]) (which gives p = (p$). The non-linear a model 
is solved by writing n = (er, tt) in terms of its longi- 
tudinal and transverse components (n • <p> = <j\<p \ and 
7r _L ip ). In the low-energy limit, the action ([3]) describes 
N — 1 non-interacting Goldstonc modes with propaga- 
tor G±(p) ~ 1/p 2 - The longitudinal propagator can be 
obtained from the constraint n 2 = a 1 + tt 2 = 1, i.e. 



G,| (r) = (a(r)a(0)) c » -^-^(r) ~ r-p^, (4) 

where (• • -) c stands for the connected part of the propa- 
gator and Gj_(r) ~ l/|r| d_2 denotes the transverse prop- 
agator in real space. Equation ([4| is obtained by using 
Wick's theorem. In Fourier space, we thus obtain 

G ll(P)-^ (P-°) (5) 

for d < 4 and a logarithmic divergence for d = 4. Con- 
trary to the predictions of Gaussian theory [Eq. @], 
the longitudinal susceptibility is not finite but diverges 
for p -> when d < 4 @, [M [H. This divergence 
is weaker than that of the transverse propagator for 
all dimensions larger than the lower critical dimension 
di, = 2. The appearance of a singularity in the longi- 
tudinal channel, driven by transverse fluctuations, is a 
general phenomenon in systems with a continuous bro- 
ken symmetry Q . The momentum scale kc below which 
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the Gaussian approximation ([2]) breaks down is exponen- 
tially small for d = 4 (and A — > 0) and of order A 1 '^ 4- ^ 
for d < 4 (sec Appendix IA 31 for the estimation of ko in 
a Bose superfluid). 

The same conclusion can be drawn from the NPRG 
analysis of the ordered phase of the linear model |T|). 
The NPRG predicts the coupling constant to scale as 
A ~ fc 4_d where k is a running momentum scale [3l| . 
This scaling follows from the flowing of the dimcnsionlcss 
coupling constant A = Afc d ~ 4 to a finite value A* for k — » 
(A ~ 1/lnfc for d = 4). The longitudinal propagator 
then diverges as Gy(p = 0) ~ ~ l/fc 4_rf and, 

identifying fc with |p| to extract the p dependence of the 
propagator, we reproduce ((5|). Thus, the divergence ([5|) 
of the longitudinal susceptibility is a consequence of the 
fixed point structure of the RG flow in the ordered phase 
of the linear model fl} . 

These considerations easily generalize to a quantum 
model with the Euclidean action 

s = \[ dT I ^{(vv) 2 +c- 2 (a rV ) 2 +V+^(v 2 ) 2 }, 

(6) 

where r is an imaginary time varying between and the 
inverse temperature (3 = 1/T, and c the velocity of the 
Goldstone mode. At zero temperature, we expect the 
(Euclidean) propagator to behave as 



G±(p,iu) 
G\\(p,iu>) 



UJ 2 + C 2 p 2 ' 



(p,w->0) (7) 



(uj 2 +c 2 p 2 )( 3 - d )/ 2 



where a; is a Matsubara frequency. (The divergence of 
C || is logarithmic in three dimensions.) The expression 
of Gii follows from with |p| replaced by (uj 2 + c 2 p 2 ) 1 / 2 
and d by the effective dimensionality d + 1 to account 
for the imaginary time dependence of the field. As in 
the classical model (fl]), it can be justified cither from 
an effective low-energy description based on the (quan- 
tum) non-linear a model or directly from the linear model 
(|6]). After analytical continuation iu> — > ui, the trans- 
verse propagator G± exhibits a pole at u> = ±c|p|. On 
the contrary, G\\ has no pole-like structure but a branch 
cut which yields a critical continuum of excitations ly- 
ing above the Goldstone mode energy uj = c|p|. This 
continuum results from the decay of a normally massive 
amplitude mode with momentum p into a pair of trans- 
verse excitations with momenta q and p q [32] . 

Interacting bosons are described by a complex field 
ip or, equivalcntly, a two-component real field (rpi, tp2)- 
In the ordered phase, the global U(l) symmetry [33[ is 
broken, giving rise to a gapless (Goldstone) phase mode 
(the Bogoliubov sound mode). Although the action of 
non-relativistic bosons differs from the relativistic-type 
action (|6j) (see Eq. (fl"3|) below), the preceding conclusions 
regarding the longitudinal propagator still hold in the su- 
perfluid phase. The reason is that the argument leading 
to ([7]) relies on the existence of a Goldstone mode with 



linear dispersion u> = c|p| rather than the precise form 
of the microscopic action. The one-particle propagator 
in the superfluid phase is usually expressed in terms of 
a "normal" self-ener gy, S n (p, iui), and an "anomalous" 
one, £ an (p,zu;) [H, H3|. In Appendix E] we show on 
general grounds that 



G\\{p,iuj) ~ - 



1 



2E an (p, iui) 



(p,w-»0). (8) 



Comparing with |[T|). we conclude that the anomalous 
self-energy 

E an (p, iW )~( W 2 + c 2 p 2 ) (3 - d)/2 (p,o;->0) (9) 

is singular at low-energy for d < 3 (the singularity is log- 
arithmic when d = 3). This singularity, which also shows 
up in the normal self-energy, is related to the infrared di- 
vergences that were encountered early on in the perturba- 
tion theory of interacting boson systems [!, H, @| ■ The 
exact result E an (0, 0) = and the asymptotic expression 
(|9|) were first obtained by NN from a field-theoretical 
(diagrammatic) approach 0, [^, |T(| ■ NN's analysis shows 
that the infrared behavior markedly differs from the pre- 
dictions of the Bogoliubov theory (E an (p, iu>) = const). 
One can estimate the momentum "Ginzburg" scale kc 
below which the Bogoliubov approximation breaks down 
from perturbation theory [l5l , l36j | (see Appendix [S| , 



(gmk h) 



)l/(3-d) 

fc^exp {-^f 

n f y gmk h 



(d<3), 
(d = 3). 



(10) 



In three dimensions, k& vanishes exponentially when the 
dimcnsionless interaction constant gmkh — > 0. In two 
dimensions, the vanishing of fc^ with gm is only linear. 

It was realized by Popov that a phase-density represen- 
tation of the boson field -0 = \fne %e leads to a theory free 
of infrared divergences [ly, [ll[ ■ Popov's theory is based 
on an hydrodynamic action S[n,9] and is valid below 
a characteristic momentum kg. Since the long-distance 
physics is governed by the Goldstone (phase) mode, a 
minimal hydrodynamic description would start from the 
phase-only action 



S =^L f dT [ d d r [ (V0) 2 +c -2 (9t0)2 ] 
2to Jo J 



(11) 



where n s is the superfluid density. This action can be 
obtained from the hydrodynamic action S[n, 9] by inte- 
grating out the density field. It is equivalent to that 
of the non-linear a model in the 0(N) model [Eq. ©]. 
Writing ip ~ y/ne 16 and expanding with respect to phase 
fluctuations (with the boson density n = const) , one finds 



G±(tt) = nGeeirr), 
G||(rr) = |G ee (rr) 2 



(12) 



for the propagator of the transverse (dip = i^fnff) and 
longitudinal (Sip = —yfn6 2 /2) fluctuations, respectively. 
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Gge is the phase propagator whose Fourier transform 
(m/n s )(p 2 + oJ 2 /c 2 ) -1 is read off from (fTTj) . fn Fourier 
space, equations (fT2|) coincide with ([7]). Thus Popov's 
approach reproduces the infrared behavior §§§ obtained 
by NN [T3| • The determination of the characteristic mo- 
mentum ko below which the hydrodynamic approach is 
valid is non-trivial in the Popov approach as it requires to 
integrate out all modes with momenta |p| > kg to obtain 
the low-energy hydrodynamic description [ll|, [l2j • Inter- 
estingly, ko coincides with the Ginzburg scale kc [IB] ■ 



II. THE AVERAGE EFFECTIVE ACTION 

We consider interacting bosons at zero temperature, 
with the action 



S 



dx 



ip*(x) [d T - (i 



V 2 \ 
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Mx)\ 4 



(13) 

(H = ks = 1 throughout the paper), where ip( x ) is 
a bosonic (complex) field, x = (r, r), and J dx = 

dr J d d r. r € [0, /?] is an imaginary time, — ► oo the 
inverse temperature, and fi denotes the chemical poten- 
tial. The interaction is assumed to be local in space and 
the model is regularized by a momentum cutoff A. We 
assume the coupling constant g to be weak (dimension- 
less coupling constant gmn 1 ~' 1 l d <C 1, with n the mean 
density) and consider a space dimension d larger than 1. 
It will often be convenient to write the boson field 



(14) 



in terms of two real fields ipi and if)2- 

To define the average effective action [3l|, we add to 
the action (Ti"3|) a source term — J dx(J*ip + c.c.) and an 
infrared regulator 



dxdx'^*{x)R k {x-x')^{x') (15) 



which suppresses fluctuations with momenta and ener- 
gies below a characteristic scale k but leaves the high- 
momenta/frequencies modes unaffected. The average ef- 
fective action 



= - In Z k [J* , J] + / dx[J* {x)<p{x) + c.c] 



(16) 



is defined as the Legendre transform of — In Z k [J* , J] , 
(Z k [J* , J] is the partition function) minus the regulator 
term AS k [<f)*, 4>]. Here cj)^(x) = (^)) J*,J is the su- 
perfluid order parameter. 

The effective action (fT5|) is the generating functional 
of the one-particle irreducible vertices. The infrared reg- 
ulator R k is chosen such that for k = A all fluctuations 
are frozen. The mean-field theory, where the effective ac- 
tion T\[4>*, (j)] reduces to the microscopic action S[4>* , 0], 



becomes exact thus reproducing the result of the Bogoli- 
ubov approximation [See Eqs. p6|) below]. On the other 
hand for k = 0, provided that R k= o vanishes, T k gives 
the effective action of the original model (fTB")) and allows 
us to obtain all physical quantities of interest. In practice 
we take the regulator 



Rkip) 



ZA.k ( 2 






LO 2 \ 


2m I* 






h k 2 c 2 ) 



(17) 



where r(Y) = (e Y — l)^ 1 and p = (p,iuu). The k- 
dependent variable Z^ k is defined below. A natu- 
ral choice for the velocity c would be the actual (fc- 
dependent) velocity of the Goldstone mode. In the weak 
coupling limit, however, the Goldstone mode velocity 
renormalizes only weakly and is well approximated by 
the fc-independent value cq = y gn/m (n is the mean 
boson density). The regulator (fT7|) differs from previ- 
ous works where R k (p) was taken frequency indepen- 
dent [H, [13, El El. The motivation for the choice (T7|) 
will appear clearly when we will discuss the BMW NPRG 
scheme. 

We are primarily interested in the effective potential 



(18) 

q> const 

(V is the volume of the system) and the two-point vertex 



T?){x,x' ;<))) = 



54>i(x)S<t>j(x') 



(19) 



<f> const 



computed in a constant, i.e. uniform and time- 
independent, field. To alleviate the notations, we now 
drop the k index. We consider = (</>i,(/>2) as a two- 
component real field [see Eq. (Tl4|) ]. U and Iy are 
strongly constrained by the global U(l) invariance of the 
microscopic action (fT3|) [33| . The effective potential U(n) 
must be invariant in this transformation and is therefore 
a function of the condensate density n = \ {4>\ + 4>2)- The 
actual (fc-dependent) condensate density no is obtained 
by minimizing the effective potential 



U'(n ) = 0. 



(20) 



Equation (|20|) defines the equilibrium state of the system. 

(2) 

On the other hand. T )■ must transform as a tensor when 
the two-dimensional vector (<pi,4>2) is rotated by an ar- 
bitrary angle a. Since one can form three tensors, 5 



-i -j 



and <pi4>j, from the two-dimensional vector U 
the most general form of the two-point vertex is 

r ij (PI 0) = Si,jT A (p; n) + <f) i (j) j T B (p; n) + e^Tcip; n) 

(21) 

in Fourier space. e,j denotes the antisymmetric tensor. 
In addition, parity and time reversal invariance implies 



rgJ(p,»«;^) = rgJ(-p I iw;^), 

r^f (p, ioj; 0) = (26 itj - l)Tf (p, -iw; <j>*), 



(22) 



5 



where <j) = (0i,02) and (jf 
(l22l) we deduce 



6 2 ). From (JST) and 



r A (p; n) = r^(-p; n) = T A (p, -iw; n), 

T B (p;n) = T B (-p;n) = r s (p, — iu>;n), (23) 

T c (p;n) = -T c {-p;ri) = -T c (p,-iw;ri). 

For p = 0, we can relate the two-point vertex to the 
derivatives of the effective potential, 

rf(p = 0; 0) = = 5i d U'(n) + (24) 



so that 



(25) 



r A (p = 0;n) = U'(n), 
T B (jp = Q-n) = U"{n), 
T c (p = 0;n) = 0. 

For /c = A, one has Tfc[0] = and therefore 



tt 1 \ 1 9 2 1 2,5/ \2 

U k =\{n) = -fm +-n = -^gn + ^{n-n ) , 



^A,k=A(p; n) = e p + g(n - n ), 
rB,fc=A(p;") = 9, 



(26) 



where n = n ,fc=A = fi/g- 

We can also relate the two-point vertex to the normal 
and anomalous self-energies that are usually introduced 
in the theory of superfluidity [H|, [H[ , 

E n (p) = + /1 - e p + f A (p) + n Q f B (p) - itcip), 

Ean(p) = n t B (p) 

(27) 

(see Appendix [A"]) , where F Q (p) = r a (p;no) (a = 
A, B, C) denotes the two-point vertex in the equilibrium 
state (n = no)- 

In the equilibrium state (n = no), the transverse 
part [37J Ta(p = 0) = U'(no) of the two-point vertex 
vanishes. This result is a consequence of the invariance of 
the effective action T[<p} in a global U(l) transformation 
and reflects the existence of a (gapless) Goldstone mode. 
When expressed in terms of the normal and anomalous 
self-energies ([27]) (with the condition Tc(p = 0) = 0), it 
reproduces the Hugenholtz-Pines theorem Q (|A3jl . 



III. DERIVATIVE EXPANSION AND 
INFRARED BEHAVIOR 



On the basis of the arguments given in the Introduc- 
tion, we expect the anomalous self-energy 

t k=0 . B {p)~{u 2 + c 2 V 2 )^/ 2 (p,w-.0) (28) 



to be singular in the low-energy limit (see Eqs. Q and 
(12710 . From 128]), we infer 



r B b = 0) = (7"(n )^fc 3 - d (fc->0) 



(29) 



Equation j29|) will be obtained in Sec.|TV]from the NPRG 
equations. In this section, we show that it is suffi- 
cient, when combined with Ward identities associated to 
Galilean and gauge invariances [f| [ID, [H[ , to obtain the 
infrared behavior of the propagators. 

The infrared regulator (TIT)) ensures that the vertices 
are regular functions of p for |p| <C k and \uj\/c <C fc, even 
when they become singular functions of (p, iuo) at k = 
(c ~ Co is the velocity of the Goldstone mode defined 
below) [39j]. In the low-energy limit |p|,|w|/c <C k, we 
can then use the derivative expansion 



r A (p) ~ Z A e p + V a uj 2 
f B (p) ~ [/"(n ) = A, 
f c(p) - Z c oj, 



(30) 



consistent with T A (p = 0) = and the symmetry prop- 
erties (f2"3"| . To obtain (f3"0")) we have expanded fs(p) 
only to leading order, dropping the next-order term 
Z B e p + V B u> 2 . Because of the singularity ([28]) . the co- 
efficients Z B and V B would diverge for k — > contrary 
to Z A , Zc and that reach finite values. The justifi- 
cation for neglecting the p dependence of the vertex T B 
comes from the fact that for d > 1, A = 0(k 3 ~ d ) is a 
very large energy scale with respect to f A , f 2 C for typical 
momentum and frequency |p|, \uj\/c ~ fc. The p depen- 
dence of T B (p) does not change the leading behavior of 
f B (p) ~ 0(fc 3_d ) which essentially acts as a large mass 
term in the propagators. 



A. Goldstone mode velocity and superfluid density 



The excitation spectrum can be obtained from the ze 

2 ma 

analytical continuation iuj — > lo + i0 + ), 



ros of the determinant of the 2x2 matrix . (p) (after 



detf W(p) = f A (p)[T A (p) + 2n f B (p)} + t 2 c {p) 
~2n f fl (0)f A (p)+r^(p) 



2n X(Z A e p + V A i 



(Z c u) 



(31) 



where the approximate equality is obtained using 
r B (p)~fc 3 - d , f A (p),f2,(p)~p 2 ,^ 2 and |p|,M/c«fc. 
Equation ([31]) agrees with the existence of a Goldstone 
mode (the Bogoliubov sound mode) with velocity 



Z A /(2m) \ 
V A + Z 2 /(2Xn )J 



1/2 



(32) 



The low-energy expansion (|30[) can also be used to de- 
fine the superfluid density n s . Suppose the phase 9(r) 
of the order parameter 4>(r) = y / 2no(cos 9(r), sin 9(r)) 
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varies slowly in space. To lowest order in V6>, the aver- 
age effective action will increase by 

p 

= no^ri 2) (p,w = O)0(-p)0(p) 
p 

ZaUq 



2m 



-p J d d r{V9f. 



(33) 



Identifying the phase stiffness with the superffuid den- 
sity [4(|, we obtain 



1l s = ZaTIq. 



(34) 



B. Symmetries and thermodynamic relations 

The two-point vertex satisfies the following relations, 

* p=0 



dp 
_d_ 



2mno ' 
1 d 2 U 



p=0 



p=0 



2n dfj, 2 
d 2 U 



(35) 



dndfi 



which follow from Ward identities associated with 
Galilean (for the first one) and local gauge (for the last 
two) invariance (see Appendix [B| . Here we consider the 
effective potential U (n, /i) as a function of the two in- 
dependent variables n and fi. The condensate density 
no = no(fj) is then defined by 



dU 
dn 



0. 



(36) 



while the mean boson density is obtained from 



'-j-U(n , p) 
dp 



dU_ 
dp 



dU_ 

dn 



driQ 
dp 



dU 
dp 



(37) 



where d/dp is a total derivative. Equation (|36[) being 
valid for any p, one deduces 



d dU 




d 2 U 




d 2 u 


dp dn 


fin 


dndp 


n 


dn 2 



dng 

no d P 



0. 



(38) 



From (f35|) and |38|) . one deduces 



V A 



Z A n = n, 
1 d 2 U 



Zo = - 



2no dp 2 
d 2 U 



(39) 



dndp 



= A 



dno 
dp 



The first of these equations states that in a Galilean 
invariant superfluid at zero temperature, the superfluid 
density is given by the full density of the fluid [f|. The 
velocity (|3"!2|) can be rewritten as 



1 



TO _d 2 U 




,a 2 u 




'dn \ 2 


dp:' 2 


n 




lit, 


\ d/j, 1 



Comparing with 

dn _ d 2 U 
du d/j, 2 

we deduce that the Goldstone mode velocity 

1/2 



d 2 U 






n 9n2 


1 

n 


\ dp J 



(40) 



^ . (41) 



m(dn/du) 



(42) 



is equal to the macroscopic sound velocity [6j . 

Since thermodynamic quantities, including the conden- 
sate "compressibility" dno /dp should be finite, we deduce 
from (|3U)l that Zc ~ k 3 ~ d vanishes in the infrared limit, 
and 



lim c = lim 

fc^O fc^o 



Z A 

2mV A 



1/2 



(43) 



Both Za = n/no and the macroscopic sound velocity c 
being finite, V A (which vanishes in the Bogoliubov ap- 
proximation) takes a non-zero value. In the infrared 
limit, the u> term of Ta(p) is therefore crucial to main- 
tain a linear spectrum and superfluidity. As discussed in 
more detail in Sec. lIV Cl the expression (|43| is a manifes- 
tation of the relativistic invariance of the effective action 
which emerges in the low-energy limit. 



C. One- particle propagator 

We are now in a position to deduce the infrared limit of 
the one-particle propagator Gy . For symmetry reasons 
(see Sec. Ell, 

Giiip; </>) = ^G u (p; n) + U 3 G tt (p; n) 



+ eijG\ t {p; n) 
for a constant field 0, where 

F^(p;n) 



(44) 



G n (p;n) = - 



D(p;n) ' 



n i \ T A (p;n) +2nT B (p; n) 
G tt (p;n) = — r , (45) 



Dip; n) 



G\tip;n) = 



Fc(p; n) 



Dip;n) 1 

and D = T 2 A + 2nT B T A + T 2 C . Equations (gSJ) follow from 
the matrix relation G _1 = -T (2) . Using A, Z c ~ k 3 ~ d > 
k 2 , we then find 



Dip) ~ 2n XV A ito 2 + c 2 p 2 ^ 



(46) 
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and 
Gn(p) 
Gu(j>) 
Git(p) 



2noA ' 
1 1 
Va uj 2 + (cp) 2 



2riQrac 2 



1 



n cj 2 + (cp) 2 ' 
mc 2 cino w 



2noAVA to 2 + (cp) 2 n d/x lu 2 + (cp) 2 

(47) 



The propagators G t t and Git have well defined limits 
when fc — > 0, while the longitudinal propagator G\\ ~ 
l/fc( 3_d ) diverges in agreement with the general discus- 
sion of the Introduction. Stricto sensu, equations 1|47|) 
hold in the limit |p|, \ui\/c <C k. We can nevertheless ob- 
tain the propagators at k = and finite (p, iw) by stop- 
ping the flow at k ~ -\/p 2 + uj 2 /c 2 (see the discussion 
in Sec. IIVC]) . Since the local gauge invariant (thermody- 
namic) quantities are not expected to flow in the infrared 
limit (Sec. IV]) . this procedure amounts to replacing n , 
n, c and dno/dfj, by their fc = values. As for the longi- 
tudinal correlation function, we reproduce the expected 
infrared singularity 



Gn(p) 



2n C[uj 2 + (cp) 2 ] 



2](3-d)/2 ' 



The constant C can be estimated by comparing ([48 
the result of Popov's hydrodynamic theory [13j, 

\ mcno / 



(48) 



with 



(49) 



From these results, we deduce the infrared behavior of 
the normal and anomalous propagators 

G n (p) = - (4>(pW(p)) 
noma 2 1 



n uj 2 + (cp) 5 
mc 2 driQ iui 



1 



GUP) 



n mc 



n dfi u 2 + (cp) 2 2 
{tp{p)tp(-p}) 

2 1 I 



Gu(p), (50) 



n or + (cp) z 2 



The leading order terms in ([50")) agree with the results of 
Gavoret and Nozieres @ . The contribution of the diverg- 
ing longitudinal correlation function was first identified 
by NN, and later in Refs. Q [H Ull • 



IV. RG EQUATIONS 

To compute approximately the effective potential U 
and the one-particle propagator, we follow the BMW 
NPRG scheme proposed in Refs. [25|, [2f| l4lf with a trun- 
cation in fields to lowest non-trivial order [42j . 



A. BMW equations 

The dependence of the effective action on k is given by 
Wetterich's equation [43[ 



9 t r[0] = iTr{i?(r( 2 )[0] +J R) x } 



(51) 



where t = ln(fc/A) and R = dtR- In Fourier space, the 
trace involves a sum over frequencies and momenta as 
well as a trace over the two components of the field <j> = 

The flow equation for the effective potential U(n) = 
(PV)~ 1 T[(f>] (with <f) = (\/2n, 0) is directly derived from 

dsr 



d t U(n) 
where 



R(p)[Gn(p;n) + G tt (p;n)}, (52) 



p J UJ 



d d p 

{2n) d 



du> 
2?r' 



(53) 



The flow equation of the condensate density is then de- 
duced from 

d t U'(n ) = d t U'\ na + U"(n )d t n = 0, (54) 

while that of A = U"(no) is obtained from 

d t X = d t U"\ no +U'"(n )d t n . (55) 

Note that the propagator G in (|52|) and below is defined 
as the inverse of — (r^ 2 ^ + R). 

Equation (|5ip leads to a flow equation for the two-point 
vertex T^ 2 ) which involves the three-point and four-point 
vertices, 

9*rgW) = 

~ 2 9 t G lll2 (q; <f>)rlf i2ii (p, -p, q, -q; 0) 



x [dtG lll2 {q;4>)]G l3ti {p + q;(f>) + {p ^ -p,i 



(56) 



where the operator dt = (dR/dt)dn acts only on the t 
dependence of the regulator R. The field 4> is assumed 
to be uniform and time independent. 

The BMW ap prox imation is based on the following 
two observations [251 ] : i) Since the function dtGij(q; <$) is 
proportional to d t R{q) 1 the integral over the loop vari- 
able q = (q, ito) in ([56]) is dominated by values of |q| and 
\oj\/c smaller than k. (Note that this argument requires 
a regulator R(q) that acts both on momentum and fre- 
quency), ii) As they are regulated in the infrared, the 
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vertices are smooth functions of momenta and fre- 
quencies [391 ]. These two properties allow one to make 
an expansion in power of q 2 /fc 2 and ui 2 /(ck) 2 , indepen- 
dently of the value of the external variable p = (p, iu> v ). 
To leading order, one simply sets q = in the three- and 
four-point vertices. We can then obtain a close equation 
for r\ (2) (p; cj>) by noting that [H[ 



rgJ(p,-p,O;0) 



r ij(p; (/>), 



_i 



rgW). 



(57) 



The flow equation for r Q (p; n) is given in Appendix [C] 
[Eqs (|C3llC5|)]. 



B. Truncated flow equations 



is more questionable in the high-frequency limit. The 

(2) 

high-frequency behavior of the two-point vertex (p) 
(and in turn of the propagator Gy(p)) follows from the 
high-frequency behavior of the propagator G(p + q) ap- 
pearing in (|5tj|) . Clearly the derivative expansion does 
not reproduce the expected high-frequency limit of the 
propagator. We shall see nevertheless that the solution 
of the flow equations does not contradict the uj v — ► 00 
limit of the propagator (Appendix |E| although the lead- 
ing corrections O(l/oj v ) and 0(l/ti; 2 ) are likely to be 
incorrect. 

These two approximations lead to the flow equations 
(see Appendix ICj) 



"tn = -h + -Itt, 



d t X = -A 2 [9J 1U1 (0) - 6Jit.it (0) + Jtt,tt( 



(61) 



and 



We simplify the BMW equations by considering two 
additional approximations. First we define the self- 
energy A Q (p; n) (a = A, B, C) by 



Ta(p; n) = e p + U'{n) + A A (p; n), 
r B ( P ;n) = U"(n) + A B (p;n), 
T c (p;n) = uj + A c (p;n). 



(58) 



It satisfies A a (p = 0;n) = 0. We then expand A a (p;n) 
about no, 



A Q (p;n) = A Q (p;n ) + (n - n )A t £> (p; n ) + ■ 



(59) 



and truncate the expansion to lowest order, i.e. we ap- 
proximate A a (p;n) by its value A Q (p) = A Q (p;no) in 
the equilibrium state. Similarly we truncate the effective 
potential U to second-order, i.e. 



U{n) = [/(n ) + ^(n-??. ) 2 , 



(60) 



where A = U"(no). For k = A, the effective action is 
given by the microscopic action S[<p] (Sec. |n}, so that 
n \k=A = fi/g, A| fc= A = g and A a \ k=A = (Bogoliubov 
approximation) . 

The second approximation is based on a derivative ex- 
pansion of the vertices and propagators. We have already 
pointed out that the integral over the internal loop vari- 
able q is dominated by small values |q|, \uj\/c < k. Fur- 
thermore, since the external variable p = (p, iui v ) acts as 

an effective low-energy cutoff, the flow of 1^ • (p; 4>) stops 
when k becomes of the order of |p| or \u>\/c. Thus all 
propagators and vertices in (|56p should be evaluated in 
the momentum and frequency range |q|, |p + q| < k and 
\ui\/c, \u> + 0J v \/c < k. In addition to the BMW approx- 
imation, we can therefore use the derivative expansion 
(f30f of the vertices in the rhs of ([56]) . This approxima- 
tion has been shown be very reliable in classical mod- 
els [!l|, 0, HE] • While we also expect a high degree of ac- 
curacy in the low-energy limit p — > 0, the approximation 



d t A A (p) = Ad t n --(/ii + 3/ tt ) 

- 2n A 2 [ Ju,tt(p) + Jttai(p) + 2 Jit.it (p)] , 

9 t A B (p) = -dtX + ^ilu-h) 
2n 

+ A 2 [-9J u>1 i(p) + J U)t t(p) 
+ Jtt.n(p) - Jtt.tt(p) + 8Ji U t(p)] , (62) 
d t A c {p) = 2n A 2 [ Jtt.it (p) - Jit.tt(p) 

- 3Jn.it (p) + 3 Jit, 11 (p)]. 

where the coefficients J a p(p) = J a fj{p',no) and I a = 
Ia{no) are defined by 

I a (n) = / d t G a (q;n), 

r „ ( 63 ) 
J a/3 (p; n) = [d t G a (q; n)]Gp(jp + q; n), 

Jq 

with a, P = 11, tt, It. The flow equations for Za, Va, and 
Zc are simply derived from 



d 

d t Z A = ——d t TA{p;n) 



d_ 

d 



dtVA = -^ d t^A{p]n) 



n—no.p—Q 



n— no,p— 



(64) 



d t Z c = -g^d t T c (p;n) 



n— no,p— 



This gives 



Z^oA 2 ^— [Ju,tt(p) + Jtt.n(p) + 2Ji t ,it(p)] „, 



5 4 Va = - 2n A 2 — [Jn.tt(p) + Jtt,n(p) + 2Jit,it(p)] p=0 , 



duo 
d - 

2n A 2 — [Jt.it (p) - Jt.tt(p) 
- 3Jii,i t (p) + 3Jit,u(p)] „. 



(65) 
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Equations (|61j) and (|65|) agree with those obtained from 
a simple truncation of the effective action T[4>] [I3|- 



C. Analytical solution in the infrared limit 

It is convenient to write the flow equations in terms of 
dimcnsionless variables 



n Q = k d Z c n 0l 

\ = k d el 1 (Z A Z c )- 1 \ 
Va = tkZ A Z7?V A 



(66) 



(see Appendix IC 3p . In the infrared limit k 0, the RG 
equations simplify, 



d t h = -(d + T)c)no, 
d t \ = (d - 2 + Vc )\ 

Q Vd+i A 

v A 



Vd+i A 2 



(67) 



Vc 



d t V A =(2 + 2r) C )V A , 
where r/c = —dt^nZc (see Appendix |P|l . We deduce 



<9 t A = (d-2)X 



and 



dtVc = (d - 2>)r]c 
For d — 3, one finds 



rib- 



Vc 



V° c 



l + rf c t' 



with rfc a constant, whereas 



(68) 



(69) 



(70) 



(71) 



rjc — > d - 3 



for d < 3. The asymptotic behavior deduced from (|70|) 
and (|71[) is summarized in Table |TJ In particular one 
finds that the coupling constant A vanishes when k — > 0, 
A ~ (lnfc) -1 for d = 3 and A ~ fc 3_d for d < 3, in agree- 
ment with the expected divergence of the longitudinal 
correlation function (Sec. IIIip . 

In the infrared limit, Zq is suppressed (Table [T]) and 
does not play any role in the leading behavior for k — > 
[Eqs. (|67|) ]. Discarding Zc, we two-point vertex r' 2 ' ex- 
hibits a space-time SO(d J r 1) (relativistic) symmetry It 
is possible to eliminate the anisotropy between time and 

~ —1/2 

space by rescaling the frequency, Q — » CjV a (the di- 
mcnsionless frequency Co is defined in Appendix IC4c|) . 
To maintain the dimensionless form of the effective ac- 
tion, one should also rescale the (dimensionless) field, 



-1/4, 



This leads to an isotropic relativistic model 





d= 3 


1< d < 3 


no 


no 


nS 


A 


(lnfc)" 1 


k e 


Zc 


(lnfc)" 1 


fc e 


Va 


VI 


V* A 


no 


(k :i Ink)' 1 


fc 2E - 3 


A 


k 


ft 1 " 6 


Va 


(kink) 2 


fc 2 - 2E 


h' 


k~ 2 


fc e " 2 


A' 


(Ink)- 1 


A'* 



TABLE I: Asymptotic behavior for fc -> (e = 3 - d). The 
stared quantities indicate nonzero fixed-point values. 



with dimcnsionless condensate density and coupling con- 
stant defined by 



V A n , 



A' 



A 



(72) 



The asymptotic behavior of n' and A' is in agreement 
with the known results of the classical (d+l)-dimensional 
0(2) model (table Hj). In particular, the dimensionless 
coupling constant A' vanishes logarithmically when d + 
1 = 4 and reaches a non-zero fixed point value A'* when 
d + 1 < 4. Using 



A = k- d (Z A e k f/ 2 V^ /2 ~X' ~ fc 3 - d A', 



(73) 



we deduce that A vanishes as k 3 ~ d when d < 3 and log- 
arithmically when d = 3. Thus, the divergence of the 
longitudinal susceptibility (which follows from the van- 
ishing of A) can be understood as a consequence of the 
low-energy behavior of the classical (d + l)-dimensional 
0(2) model. 

As explained in Appendix \D\ the infrared limit of the 
self-energies can be obtained from the derivative expan- 
sion if we stop the flow at k ~ (p 2 + uj 2 /c 2 ) 1 / 2 . This 
yields 



T A ,k=o(p) — Va^ 2 + Z A e p , 

r B ^ (p)^(Lo 2 + c 2 p 2 )^ 2 , 
rc, fc =o(p)~^ 2 + c V) (3 - d)/2 . 



(74) 



Since Z A and V A do not flow when k — > 0, they can 
be evaluated for k = and related to thermodynamic 
quantities (Sec. IIII|) . We expect the following relation 
between Tb and Tc, 



lim T _^P\ 



Z c dn 
= hm — = — — 

fe^o A a/i 



(75) 



fc=0 



10 



0.5 




— X/g 

Z c 

-- V A /V Ak=0 



2mg 




Hk G /k) 



FIG. 1: (Color online) Top panel: \/g, Zc and Va/Va.a^o vs 
ln(fcG/fc) where fee = ^/ (gm) 3 n/4iT, for n = 0.01, 2mg = 0.1 
and d — 2 [ln(fcc/fch) — —5.87]. The inset shows kc vs 2mg 
obtained from the criterion Va,Ic g = VA,k=o/2 [the green solid 
line is a fit to kc oc (2mg) 3 ^ 2 ]. Bottom panel: ln(A) and 
\n(Zc) vs \n(ka/k). A and Zc vanish as k for k <C kc (blue 
dotted lines). All figures are obtained with A = 2m = 1 and 
the regulator l|17[l. 



This relation will be confirmed numerically in Sec. [V] 
From (|T4")) and (|75|) . we reproduce the infrared limit 
(|47l48p of the propagators obtained in Sec. MI CI from 
general considerations. 



V. NUMERICAL RESULTS 

In this section we discuss the numerical solution of the 
flow equations. We consider a two-dimensional system 
in the weak coupling limit 2mg = 0.1. The actual boson 
density n = fik=o is fixed and the chemical potential /j, = 
gno.k=A is fine tuned in order to obtain n St k = ^A,fe^o,fc = 
n for k = 0. 

The flow of A, Zc and Va is shown in Fig. [TJ (The 
asymptotic behavior of various quantities as a function 
of the space dimension is summarized in table HI) In 
agreement with the discussion of Sees. [HI] and IIV1 we 
find that A, Zc ~ k are suppressed as k — » 0, while Va 
flows toward a finite value. The anomalous self-energy 
S an (p = 0) = noTsiO) = ngX therefore vanishes for k — > 




no/ n ,fc=A 
n s /n Stk=A 
c/c k=K 



Hk G /k) 
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FIG. 2: (Color online) Condensate density no, superfluid 
density n a = Zatio = n and Goldstone mode velocity c vs 
\n(kc/k). The parameters are the same as in Fig. [T] 



10.3 



10.2- 



10.1 




dfj, ~ A 



Hk G /k) 
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FIG. 3: (Color online) Condensate compressibility dno/dfi = 
Zc/X vs In(fc G /fc) [Eq. ©]. 



in agreement with the exact result Q . From Fig. [U one 
can clearly identify the momentum scale kc below which 
the Bogoliubov approximation breaks down. The inset in 
the figure shows kc obtained from the criterion Va,1cg = 
VA,fc=o/2- It is proportional to \J (gm) 3 n ~ gmkh <C kh, 



in agreement with the perturbative estimate ()10p . In 
practice, we use the definition kc = \J {gm) 3 n/Air. Note 
that the healing scale kh = \J2mgn (defined in Appendix 
|A|) keeps its mean-field (Bogoliubov) expression since the 
renormalization of the two-point vertex is very small for 
k ~ k h > k G - 

Fig.[2]shows the behavior of the thermodynamic quan- 
tities no, n s and c. Since Zam=q — 1-004, the mean boson 
density n = Zatiq is nearly equal to the condensate den- 
sity no- The condensate compressibility dng/dfj, = Z G /X 
[Eq. ([39]) ] is shown in Fig. [3l Apart from an initial vari- 
ation at the beginning of the flow {k < k^), these quan- 
tities do not vary with k. In particular, they are not 
sensitive to the Ginzburg scale k G - This result is par- 
ticularly remarkable for the Goldstone mode velocity c, 
whose expression (|32| involves the parameters A, Zc and 
Va, which all strongly vary when k ~ k G - These findings 
are a nice illustration of the fact that the divergence of 
the longitudinal susceptibility does not affect local gauge 
invariant quantities [l5j. 

In Fig. Q] we show the flow of t\a = —dtlnZA and 
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FIG. 4: (Color online) r\A and r\c vs \n{ka/k). The in- 
set shows the location of the maxima in the curves t)a (red 
crosses) and r/c (green stars) vs 2mg [the blue solid lines cor- 
respond to const x ^/2mgn oc fcjj. The parameters are the 
same as in Fig. [3] 



0.1r— 

0.08: 
0.06: 
0.04: 
0.02: 
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A| 



p\/k G ~ 1400 
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FIG. 5: (Color online) r Sifc=0 (p,iw) = A Sjfc=0 (p, iw) (sym- 
bols) and A fc _^ p 2 +U ,2/ C 2 (solid line) vs ^/o; 2 + (cp) 2 for var- 
ious values of |p|. The inset shows ln(fs i ) c= o) and a fit to 
Cy^; 2 + (cp) 2 (solid line). 



r)c = —dthiZc for k > kQ. Both 77^ and r\c exhibit a 
maximum corresponding to a slight increase of Za and 
Zc (Za then saturates to Za,u=q while Zq strongly de- 
creases when k ~ fc^). The location of these maxima is 
given by the healing scale kh (see inset in Fig. [4]) [Zrjl ]. 
The maxima of 77A and rjc become more pronounced as 
2mg increases, but remains very small in the weak cou- 
pling limit 2mg < 1. The small window around kh where 
the anomalous dimension t\a is finite is likely to be a rem- 
nant of the critical regime that exists near the Goldstonc 
regime at higher temperatures, and which is progressively 
suppressed as the temperature decreases. 



A. Self-energies 

The self-energies are obtained from the numerical so- 
lution of the flow equations ([62]) or (|C10[) . By computing 
A a (p, iuS) for TV" frequency points iui (I = 1, • • • ,N with 
typically N ~ 100), one can construct a iV-point Pade ap- 
proximant P a {p, z) which is equal to A Q (p, iu) when the 
complex frequency z coincides with one of the Matsubara 
frequencies iuj\. The retarded part of the self-ener gy i s 
then approximated by A^(p, u) = P a (p,uj + i0 + ) [471 ]. 
(All self-energies discussed in this section corresponds to 
k = 0.) 

Let us first discuss the momentum and frequency de- 
pendence of Afl(p) at k = 0. Note that Ab(p) = 
Fb(p) = £ an (p)/no s i nce ^fc=o = 0. In the following, we 
shall rather discuss Tb(p) which is the right quantity to 
consider when comparing to the Bogoliubov approxima- 
tion. Fig. [5]shows that Tb{p) is a function of to 2 + (cp) 2 , 
not only in the infrared regime \Jp 2 + to 2 /c 2 -C ka but 
also for \Jp 2 + uj 2 /c 2 > ka- Furthermore, Tb(p) is re- 
lated to the running coupling constant Afc by 

r B ( P ) - M k=v ^js- (76) 




■/(ck G ) 



FIG. 6: (Color online) Real and imaginary parts of the re- 
tarded vertex fc = o(p = 0, (J) vs Lj/(cko)- 

This confirms that Fb(p) can be (approximately) ob- 
tained from f_B,fe(p = 0) by stopping the flow at k ~ 
y/p 2 + uj 2 /c 2 . For y/p 2 + uj 2 /c 2 ^> ka, one therefore 
recovers the Bogoliubov result ^b(p) — 5, while for 
\/p 2 + oj 2 /c 2 <C ka one obtains 



T b (p) * C^ 2 + (c P ) 2 , (77) 

with C a p-independent constant. 

The Ginzburg scale ka manifests itself also in the fre- 
quency dependence of the retarded vertex r|j(p = 0, uS) 
(Fig. El). For \tu\ > ck G , the imaginary part 9[ff (0,w)] 
is very small and the real part tends to g in agree- 
ment with the Bogoliubov approximation and the exact 
high-frequency limit (Appendix [E| . But for \u\ <C cka, 
the real part is strongly suppressed and becomes of the 
same order as the imaginary part. The crossover be- 
tween the Bogoliubov and the infrared regimes can also 
be observed by varying |p| (Fig. [7]). While the Bogoli- 
ubov result Ij|(p, ui) = g is a good approximation when 
|p| ^> ka, fg(p,uj) develops a strong frequency depen- 
dence for |p| < ka- For |p| <C ka, we can use ([771) to 
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FIG. 7: (Color online) Real and imaginary parts of the re- 
tarded vertex T*g k=0 (p,u) = EaL(p,w)/no for various values 
of |p| ranging from 0.3kc up to 170/cg ~ fc/i/2 (n = 0.01 and 
2mg = 0.1). The Bogoliubov approximation corresponds to 
Ff (p,w) = s = 0.1. (K[F|] > and 9[f|] < 0.) 
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FIG. 8: (Color online) Real and imaginary parts of the 
retarded self-energies A§(p,w), A§(p,w) = f||(p,aj) and 
A§(p, u;) for |pj ~ 0.036/cg and fc = 0. The blue crosses cor- 
respond to the expressions (179[) obtained from the derivative 
expansion, while the dotted lines show the analytical result 
(I78|) obtained from the approximation (|77[) . 
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FIG. 9: (Color online) r c (p)/(wr B (p)) vs w/(c|p|) for |p| ~ 
0.036fcc and k = 0. The dashed line shows the value of 
dn /d/i|fc=o = lim fe —o Zo/A. 



obtain the low-frequency behavior (|w| <C ckc) 



fj|(p,w) ^ C v /-(^ + i0+) 2 + (cp) 2 



~ C0(c|p| - |w|)V(cp) a -w a 

-iCsgnH^I^I-clpDV^ 2- (78) 

(0(x) is the step function). The Bogoliubov result 
Tg(p,uj) = g is nevertheless reproduced for \ui\ ^> cka 
(Fig. [5]) . As shown in Fig. [5] the square-root singularity 
(|75| is also obtained from the numerical result based on 
the Pade approximant. The asymptotic result ([78]) was 
first obtained by NN within a diagrammatic approach, 
and later reproduced by Popov and Seredniakov in the 
hydrodynamic approach [ljj ■ Fig. [5J also shows the nu- 
merical results for A^(p,o>) and A^(p, u). In the in- 
frared regime |p|, \u)\/c <C he, these self-energies are very 
well approximated by their derivative expansion, 



A£(p,w) 
Ag(p,cj) 



-Vau? + (Z A - l)e p , 



(79) 



where Va =_ VA,k=o and Za = ZA,k=a- The leading cor- 
rection to A^(p,(jj) ~ iu> is given by the relation ([75)) 
between Tb(p) and rc(p), which is rather well satisfied 
when \p\,\u>\/c<k G (Fig.© 

The imaginary part of A A and the real part of Aq give 
a finite life-time to the sound mode. They arise from the 
decay of a phonon with momentum p into two phonons 
with momenta q and p q (Beliaev damping 0]). 
This damping process follows from the second contri- 
bution (proportional to f( 3 )f( 3 )) to d t T} 2) [Eq. ([56]) ]. 
Fig. [10] shows that S[A^(p,w)] and 3?[Ag(p,w)] van- 
ish for \ui\ < c|p|. The absence of damping below the 
threshold frequency ~ c|p| is due to the energy conser- 
vation u> = c|q| + c\p — q| in the decay process. While 
it appears difficult to decide from the numerical results 
whether S[Af (p, c|p|)] and 5R[Ag(p, c|p|)] (which deter- 
mine the life-time of a phonon with momentum p and 
energy c|p|) are zero or not, it is well known that for 
quasi-particles with a linear spectrum, Beliaev damp- 
ing cannot take place as there is no phase space avail- 
able [48| . Beliaev damping requires a positive curvature 
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FIG. 11: (Color online) Spectral function -i4n(p, ui) for |p| ~ 
0.036A;g and k = 0. The red solid line is the result obtained 
from the Pade approximant, while the green dashed line cor- 
responds to the analytic expression in (|8ip . The inset shows 
the ratio between A\\(p,u>) and the approximate result (|81[) 
on a larger energy scale. 



ificd from the one-loop self-energy diagrams. 



FIG. 10: (Color online) Real and imaginary parts of the re- 
tarded self-energies A^(p,cu) and A^(p,cj). For lo > c|p|, 



»[A2(p,w)] ~ a, 2 , 9f[A2(p,w)] 



5R[Ag(p, 



9[Ac(p, w)] ~ lj (blue dotted lines). The spikes are due to 
SJ[A^(p, ui)] and 5R[A^(p,oj)] being nearly zero and changing 
sign. 



of the quasi-particle dispersion, i.e. E p = c\p\ + a\p\ 3 
(a > 0). In this case, the threshold frequency, obtained 
from the condition u> = Eq + -E p - q (with p fixed), lies 
below E p . The decay of a quasi-particle into a pair of 
quasi-particles then gives a scattering rate of order |p| 3 
in a two-dimensional system [49|, [5(| . Since we use the 
derivative expansion of the vertices to compute the self- 
energies (see Sec. HV|> . the quasi-particle dispersion 
becomes linear to a very high degree of accuracy in the 
"relativistic" regime |p| <C ka- In this regime, we expect 
the curvature of the dispersion to originate in the (p, ui) 
dependence of the self-energy A^(p,w) that is not in- 
cluded in the derivative expansion. Thus a reliable com- 
putation of the Beliaev damping would require a self- 
consistent numerical solution of the flow equations. 

While the Pade approximant technique is very efficient 
to obtain A^(p, u>), as well as A^(p, ui) and A^(p, ui) in 
the infrared regime, the computation of A^(p, ui) and 
Aq(p,lj) for |p| 3> ko appears more difficult for reasons 
that we do not fully understand. (Note also that the 
use of the derivative expansion might also be a source 
of difficulties for reasons discussed in Sec. IIVBI ) In the 
limit |p| S> fee the Bogoliubov approximation is how- 
ever essentially correct and the corrections A^(p, ui) and 
A^(p, a;) provide a small broadening of the Bogoliubov 
quasi-particles (Beliaev damping) as can be directly ver- 



B. Spectral functions 

The knowledge of the retarded one-particle Green func- 
tion enables to compute the spectral functions (5lj | 



A„(p, w ) = --3[G 1 f(p,^)], 
A tt (p,u>) = —Q[G*(p,u>)], 



Ai t (p, W ) = -K[G*( P) u;)]. 

7T 



(80) 



From equations (|47|) and (|48|) . we obtain 



A u (p,u>) 
An(p,cj) 

Ai t (p,u) = i 



mcriQ 



n\p 

sgn(w) 



[6(ui - c\p\) - 6(u> + c\p\ 
uj\ -c\p\) 



2irn C ^ufl - (cp) 2 ' 
rac 2 driQ 



(81) 



2n d[i 



[5(w - c|p|) + 5(u> + c\p\ 



in the infrared regime. A t t(p,Lu) and A\ t (j),Lu) exhibit 
Dirac peaks at the sound mode energy ±c|p|. On the 
other hand, the longitudinal spectral function A\\(p, ui) 
shows a critical continuum with a singularity at the Bo- 
goliubov mode energy, in agreement with the predictions 
of the hydrodynamic approach fl3T | . The spectral func- 
tion A\\(p,<J) obtained from the Pade approximant is 
shown in Fig. 1111 The square root singularity is very 
well reproduced and extends up to lo ~ ckQ. 

From these results, we can deduce the spectral func- 
tion of the normal (U(l) invariant) Green function (see 
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Appendix [E| , 



A n (p,w) = — Q[G%(p,u>)] 

= \ Hu(p, + A tt (p, w)] - iAi t (p, w) 

~i[^u( Pj w) + Att(p,w)]. (82) 

The singularity of the longitudinal correlation function 
shows up as a continuum of excitations above the Bo- 
goliubov sound mode. The respective spectral weights 
at positive frequencies of the transverse and longitudinal 
fluctuations are given by 



S tt 
5*11 



mcriQ 

2^pT 

2 r ck o 



J— 1 (2*^ 

4™ C n \\p\J 



(83) 



for |p| <C k G . Using the estimate (|49|) of the constant C, 
we obtain the ratio 



S\\ ^ mc\p\ ( 2k G 

S t t 87TK 



(84) 



where the last result is obtained with c ~ y/gh/m. This 
ratio is extremely small in the weak coupling limit where 
mg <C 1 and |p| <C k G <C It can however become 
sizable in the intermediate coupling regime when mg ~ 1 
and k G is n °t much smaller than \J~fi. 



VI. CONCLUSION 

The BMW NPRG method provides a powerful tool to 
study interacting boson systems. In particular, it enables 
to obtain the momentum and frequency dependence of 
the correlation functions on all energy scales. Our re- 
sults reveal the crucial role of the Ginzburg scale k G in 
zero-temperature Bose superfluids. At large momenta or 
energies, |p| 3> k G or |cj|/c 3> k G , the Bogoliubov theory 
provides a good approximation to the correlation func- 
tions. For |p|, \u>\/c <C k G , the correlation functions are 
governed by a different fixed point, which corresponds 
to Popov's hydrodynamic theory. Throughout the pa- 
per, we have emphasized that interacting boson systems 
can be understood within the framework of the (quan- 
tum) 0(2) model. The infrared behavior of this model 
is characterized by singular longitudinal fluctuations in- 
duced by the coupling to transverse (phase) fluctuations, 
a phenomenon which is common to all models with a 
continuous broken symmetry Q. 



From a technical point, we have not solved the BMW 
equations in their full glory. By neglecting the field de- 
pendence of the self-energies A a (p;n) (which were ap- 
proximated by A Q (p;no)) and using the derivative ex- 
pansion, we have obtained flow equations which can be 
solved with reasonable numerical effort. Yet these equa- 
tions yield a remarkable description of the singularity of 
the self-energy induced by the divergence of the longi- 
tudinal susceptibility. Quasi-particle life-time (Beliacv 
damping) can also be obtained in principle if the flow 
equations are solved self-consistently (i.e. without rely- 
ing on the derivative expansion). 

We have restricted our analysis to the weak coupling 
limit where the two characteristic momentum scales kh 
and k G are well separated (k G -C kh <C n 1 / 4 ). The 
characteristic momentum scale n 1 ^ does not play any 
role in this limit. When the dimensionlcss coupling con- 
stant is of order unity (intermediate coupling regime), 
the three characteristic scales become of the same order: 
k G ~ kh ~ n 1 l d . The momentum range [k G ,kh] where 
the linear spectrum can be described by the Bogoliubov 
theory is then suppressed. We expect the strong coupling 
regime to be governed by a single characteristic momen- 
tum scale, namely n 1 ^. A good description of physical 
phenomena at the scale of the interparticle spacing is 
likely to require the consideration of the complete BMW 
equations (with no additional approximation) with both 
the field and (p, u>) dependence of the vertices taken into 
account. 

In one dimension, superfluidity exists without Bose 
Einstein condensation (no = 0), and our results regard- 
ing the infrared behavior of the correlation functions 
do not apply. If however, we insist on using the Bo- 
goliubov theory as a starting point, we find from the 
perturbative estimate of Sec. IA 31 a characteristic length 
k G ~ (<?m) 3 / 4 n /4 . This expression makes sense if we 
interpret n as the condensate density n ^ G at the scale 
k G . A similar characteristic scale, k s ~ (gm) 3 / 4 ??, 1 / 4 , 
has been found in Ref. [52| . In weakly interacting one- 
dimensional Bose gases, k s separates a high-momentum 
regime (|p| > k s ) where the Gross-Pitaevskii descrip- 
tion is valid, from a low-momentum regime (|p| < k s ) 
where a more elaborate description (e.g. based on the 
exact solution of the Lieb-Liniger model [HI, [54[) is re- 
quired. The description of one-dimensional superfluidity 
from the NPRG is challenging, even if the derivative ex- 
pansion yields reasonable results at weak coupling [13], 
and should be an interesting test of the BMW scheme. 
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APPENDIX A: BOGOLIUBOV'S THEORY 



self-energies 



In this Appendix, we briefly review the main results of 
Bogoliubov's theory. 



1. Beliaev's self-energies 

The action of interacting bosons is often written in 
terms of the two-component field 

*(p) = f^p ) J. &(p) = {P(p)M-p)), (Al) 

where p = (p,ico). The one-particle (connected) prop- 
agator then becomes a 2 x 2 matrix whose inverse in 
Fourier space is given by 

/ iw + \i - £ p - S n (p) -E an (p) \ 

V -Kn(p) ~ iuJ + fJ> — tp- S n (-p) J ' 

(A2) 

where S n and £ an are the normal and anomalous self- 
energies, respectively, and e p = p 2 /2m. Making use of 
(|14p and the relation G = — f ( 2 ' _1 between the prop- 
agator and the two-point vertex, one obtains equation 
(|27| if one chooses a real order parameter (ip(x)} = 
y / no- The normal and anomalous self-energies satisfy the 
Hugenholtz-Pines theorem Q 



E n (0) - S an (0) = fx, 



(A3) 



which is a consequence of the spontaneously broken 
global U(l) symmetry in the supcrfluid phase. 

Using (|A2|) . we can relate the longitudinal propagator 



= -\mp)+r{-pW{-p)+r{p)]) c (A4) 

((• • -) c denotes for the connected part of the propagator) 
to the self-energies S n and £ an . Anticipating that /i — 
T, n (p),T, an (p) ^> p 2 ,oj 2 when p,u — > (and neglecting 
terms O(p 2 ,oj 2 ), we deduce 

hm Gn{p) hm ^ _ _ _ 

-1 



lim 

p^o 2£ an (p) 

lim 



(A5) 



p^o 2n Q r B (p) ' 
where we have used the Hugenholtz-Pines theorem (|A3j) . 

2. Bogoliubov's approximation 

The Bogoliubov approximation is based on the micro- 
scopic action (|13p and a first-order computation of the 



E*(p) = 2 5 n , 



(A6) 



where the condensate density no = fx/g. This yields the 
propagators 



<£(p) = -<tf(|#*(p)) c 
<&(p) = -Wl#(-p)>c = 



3«o 



£ 2 



3«o 



(A7) 



where £7 p = [e p (e p + 2gn )] 1 / 2 is the Bogoliubov quasi- 
particle excitation energy. When |p| is larger than the 
healing momentum kh = (2gmno) 1 ' 2 , the spectrum E p ~ 
e p +gno is particle- like, whereas it becomes sound-like for 
|p| <C kh with a velocity cb = ygno/rn. In the small- 
momentum limit |p| <C kh, 



Gf^P) 



Gl(p) 



GfdP) 



j 2 + c 2 B p 2 ' 



2gn 



Ld 2 + C 2 B p 2 ' 



LU 2 + C 2 B p 2 ' 



(A8) 



Note that in the Bogoliubov approximation, the occur- 
rence of a linear spectrum is related to £ an (0) being 
nonzero. In the weak coupling limit, no is approximately 
given by the full density n, and the healing momentum 
can also be defined by kh = {2gmn) 1 ' 2 (which is the 
definition taken in Sec. |V|) . 



3. Perturbative estimate of the Ginzburg scale kc 

Let us consider the one-loop correction £ a n (p) to the 
Bogoliubov result S^(p) = gn . The leading contribu- 
tion comes from the one-loop diagram 



where the internal lines correspond to transverse fluctu- 
ations, i.e. 



E&(P) 



9 2 no 



G t t(q)G tt (p + q) 



g 2 n S d 
2 (2tt) c 



%mn,f' 2 . 



rflqllql^ 1 

2#n 



2^ \u 2 
f kh d\q 



c|q 2 



2tt) 



(A9) 



where kh is the healing momentum defined in section 
IA 21 and S<i the surface of the unit sphere in d dimen- 
sions. The infrared limit k in the integral is of order 
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(p 2 + w^/c^) 1 / 2 (with p = (p, iujy)). The one-loop cor- 
rection is divergent when d < 3. This divergence reflects 
the difficulties of diagrammatic calculations beyond the 
Bogoliubov approximation and is a manifestation of the 
diverging longitudinal susceptibility (2?| . We estimate 
the Ginzburg momentum scale fcc from the condition 

lEi^MSfLWI [sec Eq. m\- 



APPENDIX B: SYMMETRIES AND WARD 
IDENTITIES 



1. Gauge invariance 

Let us consider the microscopic action 



S 



dx 



r{x)(d T - n(x) - ^[V - xA(x)] 2 )^ 



■Mx)\ 4 



in the presence of external sources n(x) and A(x). S is 
invariant in the gauge transformation 



if>(x) -> ip(x)e ia ^\ 
iP*{x)^TP*{x)e~ la ( x \ (B2) 

fj,(x) — > /i(x) + id T a(x) 
A(x) -> A(x) + Va(x), 

where a(x) is an arbitrary real function. This implies 
that the effective action satisfies 

T[R{a)(t),n + id T a, A v + d v a] = Y[<j>, n, A v ], (B3) 

where (j> = (</>i, </> 2 ) T and 



R(a) = 



cos(a) — sin(a) 
sin(a) cos(a) 



(B4) 



is a two-dimensional rotation matrix. Differentiating 
(|B3j) with respect to ct{x), we obtain 



Est 6v 
J - r7 -^e ij (j) j {x)+td T 



54>i(x) 



Sfi(x) 



ST 



8A u {x) 



= 0. 



(B5) 

Differentiating now with respect to </>j(x 2 ) and /x(x 2 ) and 
setting </> = (y / 2n , 0), /i(x) = fi and A(x) = gives 

- v / 2raoIj2 ) (x 2 , ari ) + i<9 T1 (x 2 , x x ) 



^^ 1 fg i (x 2 ,x 1 ) = 0, 



/2nor^(xi,ar 2 ) + i9 Ti r. ( ^(a;i, x 2 ) 



(B6) 



X^ 9l 'i r ;2o( Xl -' a:2 ) = °> 



where we have introduced 



S^T 



54>i (x 2 )6(i(xi) ' 

r P), \ ^ (2) r 



(B7) 



5fj,(x2)S(j,(xi) 

and similar definitions for T?^(x2,xx) and rf 2 (x 2 , xi). 

Note that with the choice <f> = (\/2no, 0), we can identify 

to f c, and f 22 to f ^. In Fourier space, (|B6|) leads 
to the Ward identities 



Sao 


= o, 


(B8) 






(B9) 


aw 


= 0. 


(B10) 



(Bl) F rom (|B8|) . we deduce 



— r (2) f D ) 
12 (pj 



i 



p=0 



rS(p = o), 



1 d 2 U 



id (J, 



d 2 U 



n 

(BH) 



where the effective potential {/(n,/z) is considered as a 
function of both n and /j,. From (JB9J) and (jBlOjl . wc 
obtain 



— f (2) (v) 



1 5 - (2 ) 



p=0 



— r (2) ( 



; r 2 » 

and therefore 



1 



f$(0) 



p=0 



p=0 

1 d 2 U 



9 f (a 
5w 



2 r 22 V) 



^2^0 <9/i 2 

i a 2 t/ 



(B12) 



p=0 



2no cfy/ 2 



(B13) 



2. Galilean invariance 

Another Ward identity can be obtained from the 
Galilean invariance of the microscopic action. The lat- 
ter is invariant in the transformation ip{x) — > ^(x)e iq ' r , 
ip*(x) — * ^(x)e _lq r if we shift the chemical potential fi 
by q 2 /2m, which implies 



T[R(a)<f>, f i + q 2 /2m} = T[<M» 



(B14) 



where a(x) = q ■ r and the chemical potential /i is taken 
uniform and time independent. To order q 2 , equation 
(|B14|) gives 



= — ^— — — h no / dxdxT^n (% — x')a(x)a(x') 
q z ar 



.■jSVnoq^f^Cp) 
2m a^i op^ 



p=0 
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where we have set cf> = (^/2uq, 0). Since 

= -pVn 



where n is the mean boson density. 



6\l OfJL 



(see Eq. (|57| ). we finally obtain 



— f {2) (v) 



p=o 2vnn§ ' 



(B16) 



(B17) 



J 



APPENDIX C: FLOW EQUATIONS 
1. BMW equations 

In the BMW approximation, the flow equation of the two-point vertex is given by 



E 9 t G llJ2 ( 9 ;,/))rg 2n (p,-p,O,O;0) 

q.ii .12 

-■ E {r^(pA-p;0)r<^ (ci) 

q,ii—ii 

where the three- and four-point vertices in (|C1[) are obtained from the field derivatives of the two-point vertex 
[Eq. (|57|)]. From (J2TJ) and we obtain 

0) = ^[^iGiiO; n) + <5;, 2 G tt (p; n)] + e tf G?i t (p; ra), 
I« (p;<£) = ^j[rA(p;n) + (5i ) i2nr_ B (p;n)] +eyr c (p;n), 



rg( P) - P) 0;, 



|^j^,i[r^(p;n) +5i,i2nr' B (p;n)] + (<5u<5j,i + 5jj5 itl )T B (p;n) + e i j5i i iT' c (p; n)j, 



t (p, -p,0, 0;0) = -5TF{^,j^,m[r^(p;n) + ^,i2n.r^(p;n)] + {8i,i8j, m + 8j,A, m )r B (p; n) 



(C2) 



+ [<5i,l(£j,m<5*,l +^,m^,l) +^m,l(^i,;^,l + + <Si,m#j,l<*i,l] 2r7 - r s(P! n ) 

+ £i,i£j,i<^,i<5 m> i4n 2 r B (p; n) + 5j, m ey [r' (p; n) + Si tl 2nT'l;(p; n)] } 
(T' a (p\n) = d n T a (p;n), etc.) for the particular field configuration <j> = (v / 2n, 0). The flow equation (|Cip then gives 
9 t r A (p;n) = - ^I u (n)[T' A (p;n) + 2nT' A (p;n)] - ^I tt (n)[T' A ( P] n) + 2T B (p; n)} 

- 2n[Ju,tt(P5 ^)r^(p; n) 2 + J t t,u(p; n)T B (p; n) 2 + 2J lt j t (p; n)r^(p; n)T B (p; n) 

- J Uj ii(p;n)r' c (p;n) 2 - 2J 1Ut (p;ri)r' A (p;n)r' c (p;n) + 2Jn,n{p;n)T B (p;n)T' c (p;n)], (C3) 



d t T B (p; n) 



2n 



[ht{n) - h{n)]T B (p; n) - I u (n) 



5 

-T' B (p;n) +nT%(p; n) 



:ht{n)T' B (p; n) 



- J U)1 i(p; n) [X(nf + T'c(p; n) 2 ] + J n ,u(p; n) [V A (p; n) 2 + T' c (p; n) 2 ] + Jtt,n(p; n)T B (p; n) 2 

- Jtt,tt{p\n)T B (p\n) 2 + 2Ji t ^t(p;n)T B {p;n)[X(n) + T' A (p;n)] + 2J Ui i t (p; n) [X(n) -T' A (p;n)]T' c (p;n) 
+ 2Ji t .n(p; n)T B (p; n)T' c (p; n) + 2Ji t)t t(p ; n)r B (p; n)T' c (p; n), (C4) 

a t r c (p; n) = - i/ u (n) [r' c (p; n) + 2nT&(p; n)] - i/ tt (n)r^(p; ») 

- 2n| Ju,u(p; ra)X(n)r' c (p; n) + Ju,tt(p; n)V A {p\ n)T' c (p; n) - J«,it(p; ji)I\b(p; n) 2 

+ Jit,il(p; "-)r B (p; n)T' A (p; n) + Jiu t (p; n) [X(ri)T' A (p; n) - T' c (p; n) 2 ] - Ji t ,u(p; n)X(n)r B (p; 

(C5) 



where X = r^ + 2r B + 2r;,r B . The coefficients I a (n) and (p, 0), we reproduce the flow equations of the classical 
Jap(, n ',p) are defined in If we set Tc = and p = 
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0(2) model derived in Ref. jilj . 



2. Truncated flow equations 

The flow equations simplify considerably when the field 
dependence of the self-energy A a (p;n) is neglected and 
the effective potential U (n) expanded about hq as in 
(|60| . In this case the only non- vanishing field deriva- 
tive is T' A {p\n) = A while T' B (p;n) = T' c (p;n) = [see 
Eqs. (|58[)]. so that we obtain 

d t T A ( P ;n) = - -Jii(n)A--J t t(n)[A + 2r B (p;n)] 

- 2n[Ja,tt(p; n)\ 2 + Jtt,n(p; n)T%(p;n) 
+ 2Ji t) it(p; n)XT B {p; n)] , 

d t T B (p;n) = — [I tt (n) - h(n)}T B (p; n) + J lltt (p;n)X 2 

- Ju,u(p; n)[A + 2r B (p; n)] 2 

+ [Jtt,u(P; n ) ~ J tt,tt(p; »i)]r|(p; n) 
+ 4 Jh, it (p; n)r B (p; n)[A + r B (p; n)], 
d t T c (p;n) = 2n{j ttt \ t (p;n)T B (p;n) 

- Jit,tt(p; n)Ar B (p; n) 

- Jii,it(p; n)A[A + 2r s (p; n)] 

+ Jt, n(p; u )Fb(p; n)[\ + 2T B (p; n)]}. 

(C6) 

We can finally deduce the flow equations for the self- 
energy A Q (p) = A a (p;n ) from its definition (|58j). 



d t A A {p) = d t T A (p;n)\ no + r A (p;n )d t n , 

d t A B (p) = d t T B (p; n)\„ + T' B (p; n )d t n - d t X, (C7) 

d t Ac{p) = d t Tc{p;n)\ no + T' c (p;n )d t n , 

where 

F' A (p;no) = U"{n a ) = A, 

r B (p;n ) = U"'(n o ) = 0, (C8) 

T' c (p;n ) = 0. 

This leads to equations ([52^1 and (p5j) . 



3. Dimensionless flow equations 

For numerically solving the flow equations, it is useful 
to introduce the dimensionless variables as well as 
the dimensionless self-energy, 



A A (p) = (Z A e k )- l A A (p), 
A B {p) = {k- d t h Z A ZcY l A B {p), 
A C (P) = {Z A e k )- l A c (j>), 



where tk = k 2 /2m. In dimensionless form, equations 
(161162165)1 read 

3 ~ 1 ~ 
d t n = - (d + ric)n + + -I t t, 

d t ~X= (d - 2 + T] A + r] C )X 

X 2 [9Ju,u(0)-6Jit,it(0) + Jtt,tt(0)], 



2n X 2 —[j n ,tt(p) + Jtt,n(p) + 2Ji t ,it(p)] p=0) 
= - 2n X 2 — [Jtt.lt (p) - Jt,tt(p) 



- 3Ju,itO) + 3 Jit,u(p)] p=0 , 

fltVA = (2 - fj x + 27 ?c )Fa 
d 

- 2n A 2 — [J n ,tt(p) + Jtt,u(p) + 2Jit,it(p)] p=0 , 

(CIO) 

and 

d t A A (p) = ( VA - 2)A A (p) + X(I n - Itt) 

- 2n A 2 [ Ju,tt(p) + Jtt,u(p) + 2Jit,it(p)] , 

d t A B (p) = (d-2 + VA + Vc)&b(p) + ^-(At - A) 

2n 

+ A 2 [-9J 1U i(p) + Jiut(p) 
+ Jtt,u(p) - Jtt,tt(p) + 8Ji t ,tt(p)] (CH) 
+ A 2 [9 Ju,ii(0) - 6Jit,it(0) + Jtt,tt(0)] , 
«9 f A c (p) = (»M - 2)A c (p) + 2n Q X 2 [J tt ,it(p) - Jit,tt(p) 

- 3Jii.it (p) + 3Jit,n(p)] , 

where 

y = i, * = ^u, (ci2) 

and 77^ = — dth\Z Al rjc = — ftlnZc. The coefficients I Q 
and J a p(p) are defined in section [C 4 cl 

4. Coefficients I a and J a p(p) 

a. I a and J a f}{p) 

To compute the coefficients I a — I a (n ) and J a f}{p) = 
J a p(j>',no) and their derivatives with respect to p or u>, 
it is convenient to introduce 



(C13) 



(C9) 



A(p)=r A (p) + R(p), 
B{p)=A{p) + 2 m f B {p), 
C{p)=t c {p), 
D(p)=C( P ) 2 +A(p)B(p). 

With these notations, we have 



Gn = ~, G tt = ~ Gi t = §, (C14) 
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and 



where 



dtGn = -R 
dtGu = —R 
dtG\t = —R 



C 2 -A 2 
D 2 ' 
^ C 2 -B 2 

C(A + B) 
D 2 



R= -Z A e k Y( VA r + 2Yr'), 
Y 



P 2 LU 2 

H 1 

/\ L-q i\j 



(C15) 



(C16) 



with r = r(Y) and r' = dr/dY. Equations (|C14|) and 
(|C15[) can be used to compute I a and J a p{p), as well as 
cLJa/3(»U=o and <9 W 2 J a /3(p)U=o- 



This gives 

'J a 



JM = A d + 2 Z A / / dy/l 2 
de n P =o d J^J 

x { k 2 Y(r lA r + 2Yr')-? R G' a 
[VAT + (va + 4)Fr' + 2Y 2 r"]^G' a }G' p . (C24) 



The function G' a (q) can be expressed as 
1 



G'n 
GU 



D 2 
1 

D 2 
C 



(C 2 A' - A 2 B') 



(C 2 B' - A'B 2 ^ 



G[ t = -^(A'B + AB'), 



(C25) 



b. 9 £p J a /3(p)| P =o 



where A' = 9 q 2yl and B' = dq^B. Using 



Using 

G a ( P + q) = G a (q) + (p 2 + 2p • q)G' a (q) 

+ 2(p.q) 2 G^( (Z )+0(| P | 3 ), (C17) 

for p = (p, 0), we find 

J a/ 3(p) = J (d t G a ) (Op + p 2 G' p + ~P 2 q 3 G|s) + 0(| P | 4 ), 

(C18) 



and 



= 4v d 



where we use the notation 



r' — r 



(C19) 



(C20) 



(note that G a (p) is function of p 2 ). We have introduced 
v d = [2 d+1 n d / 2 T(d/2)}- 1 . U sing t he variable x = q 2 and 
integrating the last term of (|C19[) by part, we find 



d_ 

dp 



p=0 



t'tl 



d\q\\q\ d+1 (d t G' a )G' p 



The operator dt is defined by 



where 



R' = ^(r + Yr'), 
2m 

R' = -^[vat + (r)A + 4)Fr' + 2Y 2 r"] 
Zm 



(C21) 
(C22) 

(C23) 



d 


d 


d 


dR 


~ OA H 


' w 


d 


d 


d 


dR 1 


~ al 7 


f OB 7 



(C26) 



we obtain 



^G' n = 4j [C 2 {AB' + A'B + AA') - A 3 B'] , 



D 3 

^Gi u = ^[C 2 (AB> + A>B + BBi)-A>B% ^ 
d ^ C 



M G[ t = -—[ {C 2 -AB)(A' + B') 



and 



2A 2 B' - 2^'B 2 ] , 



j Q' n = -±(C 2 -A 2 



dR' 

dR' * 
_d_ 

dR 



-,G[t 



D 2 
1 

D 2 
C_ 

D 2 



(C 2 - B 2 
(A + B). 



(C28) 



Equations (|C25|) . (|C27|1 and (|C28|) are used to compute 
d €p Jap(p)\p=o- 111 the derivative expansion, we use the 
simplified expressions 

A(p) = V A uj 2 + Z A e p + R(p), 

B(p) = A{p) + 2n Q X, (C29) 

C(p) = Z c u, 



and 



A' = B' = — (1 + r + Yr'), (C30) 
2m 
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c. I a and J r a f3(p) 
We introduce dimensionless propagators, 



G tt 
Git 



Z A tk 

Gtt 



D 
B 

z A e k b ' 

Git _ G 
~Z A T k _ W 



(C31) 



where 



A = {Z A e k )- x A, 
B = (Z A e k )- 1 B, 
C = {Z A e k )- x C. 



(C32) 



The dimensionless coefficients I a and J a /3(p) are then 
defined by 



I a = k- d z c i a 



and 



2v d y d ' 2 - 1 {i lA r + 2Yr') 



,,dG a 



Or 



Z A Z C £k 

- — / d6sm d - 2 6 



x / y d,2 ~ 1 {iur + 2Yr') 



r .dG a (q)~ 



2v d 



dr 



Gp(p + q), 



y d / 2 ~\r lA r + 2Yr')^G f , 



(C34) 



(d = 3 or d = 2). y and <D arc defined in (IC12[) . To 
compute (|C34|) . we use 



9G » t { A 2 -C 2 ), 



Or 

e>G tt 
dr 

ggit 
dr 



D 2 
Y 



D 2 



(B 2 -C 2 ), 



(C35) 



£> 2 V ' 



In dimensionless form, Eq. (|C24|) becomes 



p=0 



p=0 



// 2 {(^r + 2F/)^ 

d Jy,S, 1 9r 

+ [vat + {r) A + A)Yr' + 2Y 2 r"\ 
id& 



x Y~ 



Or 



r-\-Yr' 



(C36) 



where (G' a = d y G a ) 
G'n = 
G'u = 
G'n = 



-^(C 2 A' ~A 2 B'), 

— }—(C 2 B' - B 2 A'), 
D 2 ' 

— —(A'B + AB') 
D 2 ' 



(C37) 



OG' 



Or 
dG[ t 



Or 
OG' 



dr 



r-\-Yr' 



r+Yr' 



r-\-Yr' 



2Y 
2Y 
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C 2 {AB' + A'B + AA 1 ) - A 3 B 



C 2 (AB' + A'B + BB') - A'B 3 



YC r 

^-[(C 2 -AB)(A'+B') 



2A 2 B' - 2A'B 2 



(C38) 



(C33) and 



dG[\ 
dr' 

dG' tt 
dr' 

dG^ 
dr' 



_21(c 2 -i 2 ), 

D 2 ^ ' 

^1L { C 2 -B 2 ), 
D 2 ^ ' 

™ { A + B). 
D 2 



(C39) 



We have i ntroduced A 1 = d y A and B' = d y B. In (jC36l) 
and (|C38|) . the derivative d/dr is taken with r + Yr', i.e. 
R' = dq2R. fixed. If A, B and C are evaluated within 
the derivation expansion [Eqs. (|C29IC30[) ], 



and 



A = V A u 2 +y + Yr, 
B = A + 2n A, 
C = u>, 



A' = B' = 1 + r + Yr'. 



(C40) 



(C41) 



APPENDIX D: SOLUTION OF THE FLOW 
EQUATIONS IN THE INFRARED LIMIT 

In this appendix, we consider the regulator ((TT] 
with (El 



We also take 



r(F) = i_^(l-F). 



Z A 
2mV A 



(Dl) 



(D2) 
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and note that Co is fc-indcpendent in the infrared limit 
k <C ka and equal to the the Goldstone mode velocity c 
(Sees. HJ and©. For Y < 1, one then has 

A = Y + Yr(Y) = 1, 

£? = 1 + 2n A, (D3) 

£> = 1 + 2n A + o; 2 . 

We also observe that the condition Y < 1 implies 



p| < ^c fc ~ fc 2 ~ d , 



(D4) 



where we have anticipated that Zc ~ fc 3 d for d < 3. On 
the other hand, 

n A = (^ ACfc ) _1 noA ~ fc 1 "'. (D5) 

We can therefore neglect <I> 2 with respect to J3, and 

D~B~ 2h X (D6) 

becomes frequency independent. For d = 3, \uj\ < 
l/(fclnfc) and n A - l/(fc 2 lnfc), so that i[D6]) holds. 

We are now in a position to compute the infrared limit 
of the coefficients I a and J a| g. Since tja — > 0, we have 



(D7) 



where we have used the fact that the condensate density 
no flows to a finite value when fc — > (so that the flow of 
no is determined by the purely dimensional contribution). 

With a similar reasoning, we find 
d t \~ (d-2 + 7 lc )\~ A 2 J tt ,tt(0) 



Vc 



(d-2 + r lc )\ + 8 

~ d ~ 
-2h \ 2 — Jtt,it(p) 



^rf+i A 2 
d + 1 V \ /2 ' 



d t V A ~ (2 + 27 7c7 )Va. 



(D12) 



L> 2 



All the integrals involved in the derivation of 
(1D11ID12[> are d + 1-dimensional integrals of the type 
(|D9p . This is a direct manifestation of the rclativis- 
tic invariance which emerges in the low-energy limit 
(Sec. IIV C[) . To compute the infrared limit of the flow 
equations satisfied by the self-energies, wc need to com- 
pute the coefficients J a /3{p) for finite p. The external 
variable p acts as a low-energy cutoff, so that J a p (p) can 
be obtained from J a /3(p = 0) with fc ~ (p 2 + oj 2 /c 2 ) 1 / 2 
(this choice satisfies the relativistic invariance). 



APPENDIX E: HIGH-FREQUENCY LIMIT OF 
THE TWO-POINT VERTEX 



Since \u)\,A -C B, we can neglect la with respect to I tt 
and approximate 

J y.Cj 'J 

= ±v d f y^Oil - Y). (D8) 
For any function f(Y), 

/>00 />00 1~ 

JO J-oo 27T 

= VX 1/2 v d +i r dY Y^-WfiY), (D9) 



so that we finally obtain 



Vd+l f>-l/2 

d + l 



I tt * 8^V A 



(D10) 



and 



d t n ^ ~(d + r) C )n + -I tt 

/ j I x~ , , Vd+l fV- 1/2 

~ -(d + r)c)n +4-—V A 
~ -(d + r)c)no, 



(Dll) 



The normal and anomalous propagators G n and G a 
defined in Appendix [A] can be written as 



Gnip) = 
G an (p) = I 



ILU — Lj' 

, A an (p,o/) 



(El) 



when iw ^ 0. The spectral functions A n and A an are 
defined by 



A n ( P ,t) = i-([vi(p J t) ) Vit(p ) o)]), 

4*>(P,*) = i([Vi(p,i),Vi(-p,0)]), 



(E2) 



where V>(p, i) and ?M(p,i) are the boson operators in 
the Heisenberg picture. From the spectral representation 
(|Eip . we obtain the high- frequency expansion 



G n (p) = - + -% + 0( W - 3 ), 

G an (p) = -^ + 0(Lu- 3 ), 
(io;) z 



(E3) 
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where 

/oo 
duicoA n (p, u>) 
-oo 

= 2ir [id t A n (p,t)] t=0 = H], ft (p)]}, 

/OO 
duju;A an (p,u;) 
-oo 

= 2tt [id t A m (jp,t)] t=0 = ([[4>(p),H],i>(-p)]), 



(E4) 



and _ff is the quantum Hamiltonian corresponding to the 
action (fT3|) , To obtain (|E4|) . we have used the equa- 
tions of motion of the operators ip(p,t) and ft(p,t). A 
straightforward calculation gives 



X r 



fi + 2gn, 



(E5) 



where n = (ft (r)ip(r)} is the mean boson density. In- 
verting (|A2|) and considering the high-frequency limit, 
we obtain 



X p = - lim [fx-e p - E„(p)] , 
Fp = - lim E a „(p), 

CJ — s-OO 

lim E„(p) = 2#n, 
lim E an (p) = 3 (V3(r) 2 ). 



(E6) 



(E7) 



From d27J) and ([581), we finally deduce 



= lim A A (p) = -M- .g(V^r) 2 ) + 2 5 n, 



f f = lim f B (p) 



(V>(r) 



(E8) 



A; 



lim A c (p) = 0. 



From these limiting values, we can obtain the "pairing" 
amplitude (^(r) 2 ) and the mean boson density n. In the 
weak-coupling limit, n ~ no ~ (ip(r) 2 ) and no ~ /i/p 
does not differ much from the Bogoliubov result, so that 
we expect A% < gn, f f ~ g and Ag 5 = 0. 

Since lim^^oo J a p{p) = 0, the flow equations (|6II62D 
yield 



dtfS 



A (In 
A 
2n 
0. 



Iu - hi) 



(E9) 



These equations arc not exact as they involve A ~ 
r^ 4 )(0, 0, 0, 0) rather than the high-frequency limit 
lim^^oo T (4) (p, —p, q, -q) (with p = (p, «£*;„)) of the 
four-point vertex. Nevertheless, the numerical results of 
Sec. [V] arc in good agreement with the asymptotic values 
(|E8|) . Note that contrary to ([E9]) . the BMW equations 
would be correct in the high-frequency limit. 
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